Variational wave functions for homogenous Bose systems 
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We study variational wave functions of the product form, factorizing according to the wave vectors 
k, for the ground state of a system of bosons interacting via positive pair interactions with a positive 
Fourier transform. Our trial functions are members of different orthonormal bases in Fock space. 
Each basis contains a quasiparticle vacuum state and states with an arbitrary finite number of 
quasiparticles. One of the bases is that of Valatin and Butler (VB), introduced fifty years ago and 
parametrized by an infinite set of variables determining Bogoliubov's canonical transformation for 
each k. In another case, inspired by Nozieres and Saint James the canonical transformation for 
k = is replaced by a shift in the creation/annihilation operators. For the VB basis we prove 
that the lowest energy is obtained in a state with ~ V volume quasiparticles in the zero mode. The 
number of k = physical particles is of the order of the volume and its fiuctuation is anomalously 
large, resulting in an excess energy. The same fluctuation is normal in the second type of optimized 
bases, the minimum energy is smaller and is attained in a vacuum state. Associated quasiparticle 
theories and questions about the gap in their spectrum are also discussed. 

PACS: 03.75.Hh, 05.30.Jp 



I. INTRODUCTION 

Since now more than half a century the understand- 
ing of Bose-Einstein condensation (BEC) in interacting 
Bose gases has been representing a major challenge for 
theoretical physics. In spite of a considerable progress, 
the decisive step in rigorous theory — simply, the proof 
that BEC indeed occurs — still could not be made. Ap- 
proximate or effective theories that take BEC for granted 
have also been struggling with some emblematic prob- 
lems as, for example, the question of the gap appearing 
almost inevitably in them. The search for gapless approx- 
imate theories has a long history that will not be evoked 
here; instead, we refer the reader to the review of Za- 
grebnov and Bru [1] and to the introduction of Yukalov 
and Kleinert's paper [2]. The most successful effective 
theory is due to Bogoliubov [3], [4], [5]. Based on the 
hypothesis that BEC takes place in the ground state, 
Bogoliubov introduced some approximations in order to 
turn the problem into soluble. He dropped part of the 
Hamiltonian (Bogoliubov truncation, BT), and replaced 
Uq and ao by volume (V) dependent complex numbers 
(Bogoliubov substitution, BS). The upshot of the ap- 
proximations is a Hamiltonian that can be diagonalized 
by Bogoliubov's canonical transformation. In the subse- 
quent physics literature BT and BS appear either sepa- 
rately or coupled together, with an embarrassing variety 
of truncation schemes. In the past and also more recently 
much effort has been made to justify these approxima- 
tions. Their status from the point of view of mathemat- 
ical rigor and also of physical significance proved to be 
quite different. 

In a fundamental paper [6] Ginibre studied BS alone 
and proved that if the pressure is maximized with respect 
to the complex number substituting ao, the approximate 



pressure converges to the true one in the thermodynamic 
limit. This still did not imply that BS was indeed "ex- 
act" , because concerning the states, it represents a dras- 
tic reduction of Fock space: only states of the form zip are 
considered, where z is a complex number, z = e^"""^"" 
and tp is an element of the Fock-subspace F' built on the 
plane waves {e'^^ '' /VV}]i^o. Justification of BS means 
therefore to prove that along with the pressure the av- 
erages of the physically most relevant operators are also 
correctly reproduced by it. A step in this direction was 
made in two recent papers [7], [8]. These showed that if 
flo — > z{V) is a complex number substitution that makes 
the pressure of the BS model converge to the true pres- 
sure in the thermodynamic limit, then |2;(y)p/T^ and the 
exact {aQao)/V have the same limits as V tends to in- 
finity. The result implies that BS predicts BEC, when it 
occurs, with the right density of condensate. 

No similar rigorous justification exists for any kind of 
truncation of the Hamiltonian. While BS alone is quanti- 
tatively correct in the sense described above, BT cannot 
be hoped to be quantitatively correct except for the cases 
of extremely low densities and weak interactions. Even 
qualitative correctness can partly be jeopardized, as it 
happens for instance in relation with the gap. Trunca- 
tion of the Hamiltonian changes the ground state energy 
in the order of the volume and introduces unphysical fea- 
tures as the non-locality of the interaction. By this, it 
causes the breakdown of the local gauge invariance of 
the field operators in coordinate space, and changes the 
global U{1) symmetry into local one in k-space. There- 
fore, truncated Hamiltonians should be considered as sep- 
arate models that can describe certain but not all ex- 
pected properties of the physical system. 

There exists a particular class of truncated Hamiltoni- 
ans in which truncation is variational in the sense that 
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the cigcnvahics are upper bounds to the lowest eigen- 
value of the physical Hamiltonian in every eigensubspace 
of the total momentum. Such a family of Hamiltonians 
appeared fifty years ago in a paper by Valatin and But- 
ler [9]. The Valatin-Butler (VB) Hamiltonians depend on 
the parameters of Bogoliubov's canonical transformation 
and are diagonal in the quasiparticle creation and anni- 
hilation operators obtained through this transformation. 
The corresponding quasiparticle vacuum state contains 
pairs of physical particles of opposite momenta. The ex- 
citation of such pairs was thought already by Bogoliubov 
to be the dominant effect of the interaction on the ground 
state of a Bose system. Valatin and Butler therefore ex- 
pected to obtain a good quasiparticle theory by mini- 
mization of the energy of the vacuum state over the set 
of parameters of the Bogoliubov transformation and by 
retaining the diagonal Hamiltonian defined by the mini- 
mizing parameters. They failed in this program, but their 
paper is remarkable and to locate the source of their fail- 
ure even fifty years later seems to us relevant. Our con- 
clusion, based on the analysis of Section III and valid for 
positive pair interactions with a positive Fourier trans- 
form, can be summarized in two points. First, within the 
family of eigenstates of the VB Hamiltonians the state 
of lowest energy is not a vacuum state but a state with 
~ \/V quasiparticles of zero momentum. Second, Bogoli- 
ubov's transformation in the zero mode leads to anoma- 
lous fluctuations of the number of condensed physical 
particles, increasing the lowest attainable energy. 

The two problems of the VB approach arc related 
and both disappear if Bogoliubov's transformation in 
the zero mode is dropped and replaced by Bogoliubov's 
substitution. This was first done by Nozieres and Saint 
James [10]. In Section IV we analyze in detail and some- 
what extend the Nozieres- Saint James (NStJ) approach. 
First, BS is replaced by a shift, a^ — z, according to Araki 
and Woods [11]. This permits us to keep the zero mode 
in consideration and to allow deviations from a coherent 
state. We then introduce variationally truncated Hamil- 
tonians defined now on the full Fock space and depend- 
ing on the parameters of Bogoliubov's transformation for 
k 7^ and on z for k = 0, consider the family of their 
eigenstates, and prove that the lowest energy is attained 
on a vacuum state with an energy density lower than the 
VB minimum. Finally, we add a Bogoliubov transfor- 
mation also in the zero mode and minimize the energy 
of the vacuum state with respect to it. While the gain 
in energy is only of the order of unity, by this way the 
gap of the associated quasiparticle theory, discussed in a 
moment, can be made continuous at k = 0. 

Although we are interested only in minimizing the en- 
ergy of our trial functions, because of their appearance 
as cigcnstatcs of quasiparticle Hamiltonians we compute 
the corresponding gaps and comment on them at sev- 
eral places. In the spectrum of a physical Hamiltonian 
with a translation invariant and not too slowly decaying 
interaction there is no gap above the ground state (in 
fact, there is a gap that disappears as V tends to infin- 



ity). This was first shown by Hugcnholtz and Pines [12] 
for k = 0, then by Bogoliubov [5] for k ^ and some- 
what later a rigorous proof, related with the breakdown 
of Galilei invariance and not with BEG, was given by 
Swieca [13]. On the other hand, truncated models most 
often produce a gap to excitations. Sometimes the gap 
disappears for k = but stays for k — > (see [14] and 
Section HI), sometimes it stays for both and the two gaps 
can be different (Sections III.B and IV. B) or can coincide 
(Section IV. G). It is because of this gap that the work 
of Girardeau and Arnowitt [15] was considered later to 
be inconsistent [12], although they minimized correctly 
the energy of their trial state. (Working with a fixed 
number of particles, they avoided the problem of anoma- 
lous fluctuations.) However, a gap does not mean the 
internal inconsistency of an approximate theory. The 
quasiparticle theories emerging as a by-product from the 
present variational treatment are, in fact, consistent by 
themselves. They can be gapful without contradicting 
Goldstonc's theorem applied to them because there is no 
BEG of quasiparticles in their ground state. 

II. PRELIMINARIES 

We are going to perform a thorough study of varia- 
tional wave functions of the product form, factorizing 
according to the wave vectors, for the ground state of 
the Hamiltonian 

H = ^EO<^)Nk - f^N + w(q)ak+q</_qak'ak 

q,k,k' 

(2.1) 

Here aj^ and ak create, respectively, annihilate a bo- 
son in the single-particle state (l/VV) exp{ik- r}, where 
k — (27r/L)(ni, n2, ns) with m integers and = V, 
£(k) = h?k'^/2m, ji and w(q) = w(— q) are real numbers, 
iVk = ajuk and N = '}2 -^k- H is defined with periodic 
boundary conditions on a cube of side L; v is the Fourier 
transform of an integrable pair potential yj, 

v{\i) = j ip{v)e-'^ '' dr, (2.2) 

and most of the time both v and ip will be supposed to be 
nonnegative. A typical example of a positive interaction 
with a positive Fourier transform is a Gaussian interac- 
tion. More examples can be obtained by convoluting a 
real nonnegative function '0(x) with -(/;(— x): if 

(^(x) = j ^(y)V(y + x)dy, (2.3) 

then both ip and v = (p = are nonnegative. Note 
that the chemical potential is included in the definition 
of H, so the ground state of H is understood in Fock 
space. For <p > the interesting region of the chemical 
potential is /i > 0; for < the ground state is the 
vacuum. 
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A. Pairing product states 



with 



Pairing product states are wave functions of the general 
form $ = $|0) where |0) is the physical vacuum state and 

$ = ( E Po,m{a*or] n ( E P^,m{ata\r] 

\m=0 / {k,-k}5^{0} \m=0 / 

(2.4) 

A variational ansatz of this kind suggests a division of H 
in three parts, 



H = Ho + H^o + 



R 



(2.5) 



where Hq collects all terms containing A'^Oi -^o; = 
and Aq, H^q collects all terms containing A^k, = 
flkfl-k and only for k 7^ 0, and Hu is the sum of the 
remaining terms of H not contributing to the energy of 
Hr^) = 0. Thus, 

Ho = -/z7Vo + ^(iV2-iVo) + ^E[^(0) + ^(k)]iVk 



+ ^E^(i^)[^s^k+^k^o] 



(2.6) 



kjiO 



and 



H^o = E[^(k) - M]iVk + ^ E(^k - iVk) 



k^tO 
1 

+ 2-y 



k^^O 



^ E KO)+«(k-k')]iVkA^k' 

OT^kT^k'T^O 

-^E E v{k-k')At,A^. (2.7) 



kT^O k'T^O.ibk 

For an operator B let (B) = ($, B$). In (ffo) and (H^o) 
the mean value of A^/Ak factorizes. This makes it possi- 
ble to choose (Ak) independent complex variables: pre- 
scribing any set {^k} of complex numbers, $ can be cho- 
sen so that (Ak) = zu for all k. The Valatin-Butler vac- 
luim state with complex gik serves below as an example. 
Define 

^1 = ^EK0) + t;(k)](7Vk) 



= ^^v{k)Re{Ay,) . 



(2.8) 



kT^O 



Choose (Aq) to be real and nonnegative. When mini- 
mizing (H) , the sign of U2 has no influence on the value 
of {H^o) but i>'2 < is necessary to minimize {Hq). We 
then find 

(^^o)= - + H + (iVo) + ^(A^o)^ 

+ ^D'iNo) + W2\[{No)-{Ao)] (2.9) 



D\No) = {N^) - {Nof . 
Observe now that 

(TVo) = ||ao1>||' , {No) + l = \\a*n' 
and therefore by Schwarz inequality 



(2.10) 

(2.11) 



{Ao) < {No)^/l + l/{No) < {No) + \ . (2.12) 

Because D'^{Nq) > 0, we obtain the following lower 
bound on {Hq): 

(Ho) + H + f^) (iVo)+f) (iVo)2-M , 

(2.13) 

The minimum in (A^o) of this bound is either — |f2|/2 if 
11-1/1 + 11/21 < or 

(Ho) > -^^{n-t^i + W2\r-l{n-i^i+2W2\)+0{V-^) 

(2.14) 



obtained with 

(No) 



+ + ^ (2.15) 



if — 1^1 -h \v2\ > 0. A particular case of (2.4) is when $ is 
an eigenstate of A^o, so among all /3k, m only (3o,m is non- 
vanishing for a single value of m. Then (A^k) = (^k) = 
for kj^O, {Ao) = D^{No) = and 



(H) = {Ho) 



2V 



m+^m^ (2.16) 



leaving m as a single variational parameter. Minimiza- 
tion then yields 



(2.17) 



obtained with m = fiV/v{0) + 1/2. This is the ground 
state energy of the mean field and also of the so-called 
perturbed mean field or full diagonal models [14, 16, 17]. 
It will serve as a reference value, the parameters will have 
to be chosen so as to obtain an energy lower than Emfi 
by getting the closest possible to the lower bound (2.14) 
with the largest possible value of 11/2! — z^i without loosing 
too much energy in (i?-,o)- 



B. Coherent product states 

The ansatz (2.4) will be partly restricted, partly ex- 
tended. We introduce two sets of quasiparticle annihila- 
tion and creation operators. 



flk = Ok — 2k, 



Zk 



(2.18) 
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and 



bt = 



1 



(S^-5l^S_k) (2.19) 



where Z],, and £fk are complex numbers, |(7k| < 1 and 
gk = 5-k- The coherent state 



Z = exp 



^(^kOk - ^^kCtk) 



|0) 



(2.20) 



is the common vacuum of each ak, a\^Z = 0. It is then 
straightforward to veriiy that 



$o = n(l-|5k|')^exp 



Skflkalk 



(2.21) 



is the vacuum state of each 6k, i-e. 6k^o = 0. On $o one 
can build up an orthonormal basis in the Fock space, 



(2.22) 



by letting {rrik} run over terminating sequences of non- 
negative integers. Unless 2;^ = for all k 7^ 0, Z and 
thus ^{mk} ^-re not eigenstates of the total momentum. 
Choosing l-J-k] = l-^k] one has at least 

(E^^^k) ^{m,}) = ^mkk . (2.23) 



Given a basis (2.22), it is possible to define a quasipar- 
ticle theory of the interacting Bose gas via the Hamilto- 
nian 

}>(^{m.}| (2.24) 

{mk>0} 



where 



{"Ik} 



($ 



{nik} 



(2.25) 



Such a theory has the merit of being variational to the 
ground state of H. If 2k = for k 7^ and, hence, 
Hqp commiitcs with the total momentum operator, the 
quasiparticlc theory is variational also in each momentum 

eigenspace. 



E{mk} > inf spec 1111^^^} 



(2.26) 



where H^^^} projects onto the subspace of states of mo- 
mentum Ti^mkk. Those concerned with effective the- 
ories of the interacting Bose gas may not optimize Hqp 
by minimizing its ground state energy, but try instead to 
make the gap above the ground state disappear. In this 
paper we insist on an unbiased minimization of the en- 
ergy, yielding what we call the quasiparticle ground state 



(QPGS) of H. At several places we shall comment on 
the gap in the spectrum of Hqp. 

If 2;k = for all k ^ 0, the dressed vacuum state 
(2.21) and states (2.22) with m_k = rny^ are also pair- 
ing. If Zk 7^ for some k 7^ 0, ${mk} provides a non- 
vanishing average of Hn and, thus, is not of the type 
(2.4). Such product states can have lower energies but 
not necessarily in the order of the volume. One spe- 
cial case of them, when all gk = 0, can immediately be 
excluded. Indeed, a simple calculation shows that when- 
ever t; > 0, m.m{Z,HZ) = -ij?V/2v{0) > Emf, obtained 
with Zk = ■\/ /iy/u(0)(5k.o- With the benefit of hindsight 
one can state that for k = the shift (2.18) is more im- 
portant than the Bogoliubov transformation (2.19) while 
for k 7^ the opposite holds true. A full variational treat- 
ment, extending to the sets {^k}, {^k} and {mk} seems 
at any rate too complicated. We will therefore focus on 
the case = Z(5k,o- The variational parameters are then 
z and gk (including 170), and all basis functions (2.22) will 
be considered. Two important special cases are z = 0, 
the Valatin-Butler (VB) scheme and go = 0, the Nozieres- 
Saint James-Araki- Woods (NStJAW) scheme. 



III. THE VALATIN-BUTLER SCHEME 

Valatin and Butler [9] , inspired by the then fresh BCS 
paper [18], proposed a quasiparticle theory for bosons in 
the lines described above, with Ck = ak (all Zk = 0) 
and (/k real. The dressed vacuum (2.21) is then pair- 
ing also for k = 0. In choosing gu real, Valatin and 
Butler followed Bogoliubov [4]. In principle, this choice 
restricts generality since the most general gauge trans- 
formation that leaves H invariant for 
all t G and a G K. is not enough to eliminate an arbi- 
trary k-dcpendcnt complex unit factor of g^. However, 
the ground state of H is unique and one may suppose 
that its QPGS is also unique. Then ^k can indeed be 
chosen to be real. 

Valatin and Butler made some ad hoc assumptions 
which eventually falsified the minimization of the en- 
ergy. Their condition on the interaction was vague; they 
thought it to be repulsive but, most important, such that 
the QPGS is the vacuum for some {gk}. This assump- 
tion was inconsistent with the variational determination 
of 9o. so they chose it to make the gap to k = excita- 
tions vanish. However, because 



($o,A^o$o) 



9'o 



9l 



(3.1) 



it was in fact the condensate density that they adjusted to 
make the gap disappear. This was a strongly contestable 
step for at least two reasons. First, the gap in question 
is in the spectrum of Hqp, therefore the breakdown of 
a continuous symmetry of Hqp could only be in conflict 
with it. If the groimd state of Hqp is the vacuum state, 
a gap above it does not contradict Goldstone's theorem 
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applied to -ffqp. One should be worried only if Hqp had 
a ground state with a macroscopic expectation value of 
6q6o and still there would be a gap above the ground 
state. Second, the density of the condensate of physical 
particles can be varied to minimize the energy, but not 
for other purposes. The only free physical parameter is 
the chemical potential If the Af-particle interaction 
energy is superstable (meaning that in any finite volume 
V it eventually increases with N not slower than N'^ /V), 
the full particle density is a strictly increasing function 
of /i (at zero temperature for > 0, if the interaction 
is nonnegative). Therefore, instead of n one may freely 



prescribe the full density, but not that of the condensate. 

The explicit form of the quasiparticle Hamiltonian 
(2.24) can be found e.g. by rewriting H in terms of the 
operators 6^ and 6^ and keeping only those terms com- 
muting with each Mk = bl^bu- Let 



2 

= ($0, iVk$o), Xk = = ($0, ^k$o) 



(3.2) 



so that Xk = + !)• A lengthy but straightforward 

calculation yields 



J 



Hqp =wo + ^ ekMk + ^ Wkk(M^ - Mk) + ^ Wkk'MkMk' 

k k kjik' 



(3.3) 



where 



Wo 



^(e(k) - M)/ik + ^ f^/ik) + oTF E^(k - k')(xkXk' + /ik/ik') 
k \ k / k,k' 



(3.4) 



ek = (1 + 2/ik) ^e(k) -1^+^^ HO) + v{k - k')] h^^ + 2xk^ " 



(3.5) 



wkk = + 2/ik)^ + Ar^Xk 



2V 



V 



(3.6) 



«^kk' = ^[f(0)+t.(k-k')] + ^[2f(0)+f(k+k')+f(k-k')](/ik+/ik'+2/ik/ik') + ^b(k+k')+t'(k-k')]xkXk' • (3.7) 



r 



It is seen that e_k = ek and w_k,-k' = ^kk' = Wk'k 
for all k, k'. From Eq. (3.3) the energies (2.25) can be 
obtained by replacing the operators Mk by the integers 
mk. Alternatively, one can take the average with ^{mk} 
of Eqs. (2.6) and (2.7) and use 

(iVk) = /ik + (1 + /ik)TOk + /ikTO-k, (3.8) 

(A) =Xk(l + mk + m_k), (3.9) 

= (A^kA^-k)-(iVk)(Af-k) 

= Xk(l+<5k,o + mk + m_k + 2mkTO_k). 

(3.10) 

The expressions for wq and Ck were already given in [9]. 
A particular member of the family of operators (3.3), 

obtained with all gk = 0, is the perturbed mean field 
or full diagonal (FD) Hamiltonian which has long been 



serving to test different ideas concerning the interacting 
Bose gas [14], [16], [17]. It will prove to be useful here as 
well. In a more convenient form it reads 

+ ^ E^(k-k')A^kAfk' . (3.11) 

k^^k' 

It is clear that the overall minimum of i?{mk}({.9k}) can- 
not be higher than the ground state energy of -ffpD- Sup- 
pose that v{li) > and is continuous at with v{0) > 0. 
Then for < the ground state of Hpu is the physical 
vacuum with zero energy. If > 0, it is (mo!)~ 5 (a*)™o |0) 
with 

mo = I^V/v{0) + ^ = mv- (3.12) 
(We may suppose that my is an integer.) This gives 

Emo=my=-^^(l + ^^ =EmF (3.13) 
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for the ground state energy of -ffpD, cf. Eq. (2.17). What 
about the gap to excitations? If /i < 0, the gap is at 
least to any excitation of the vacuum. If ^ > 0, 
because toq = my was chosen to minimize the quadratic 
polynomial — /umo + [v{0)/2V]{mQ — mo), we find 



Er, 



^(0) 2 

n 

2V 



(3.14) 



which for n = o{VV) tends to zero as V goes to infinity. 
Thus, the gap to k = excitations vanishes, and this is 
consistent with the breakdown of the gauge invariance 
during BEC in this model system [14]. On the other 
hand, for k ^ 

v(0) + v(k) vCk) , 

= e(k) - M + ^ ' ' mv > li^^ . (3.15) 



(Above it is understood that m\^i = unless otherwise 
stated.) So there remains a partial gap to k ^ exci- 
tations in the hmit of infinite volume, which tends to fj, 
with k going to zero. 



A. The quasiparticle ground state in the VB 
scheme 

Proposition 1. Let the pair interaction be nonnegative, 
> 0. Then tlic quasiparticle ground state is not a 

dressed vacuum state. If also v > 0, the energy difference 

between any vacuum state and the QPGS is of the order 

of the vohime. 

Proof. Note that {<^q,H<^q) = wq. > implies that 

V is positive definite, that is, 



k,k' 



(k-k')aa' >0 



for any choice of ^k- Therefore, even without supposing 
u(k) > 0, 



Wo 



X] £k/ik + ^ X] ^('^ ~ k')(xkXk' + /ik/ik') 



k,k' 



MO) ^ ^ 



If V > 0, the lower bound is slightly above the ground 
state energy (3.13) of the FD model. Otherwise, using 
(my!)~ 2 (a*)™!' |0) as a variational wave fmiction, one 
can see that the ground state energy of i/pD is bounded 
above by the value (3.13), so is certainly smaller than 
wq. This proves that the QPGS is not a vacuum state. 
Assume now <p > and w > 0. A lower bound on wq is 
obtained by dropping 



«(k - k')xkXk' > 



from the right-hand side of Eq. (3.4). Let = 
If fc = |k| > then £(k) - > 0. The 
other terms for or fc' > k^ being also positive, the lower 
bound on wq can further be decreased by restricting the 
summations to fc, fc' < fc^. Dropping also ^e(k)ftk > 0, 
we find 

Wo > -iJ-x + ^^^x^ + ^ "^C^ ~ k')/ik/ik' 



2V 



1 



2V 

fc ) fc' fcyj, 

v{Q) + min ^(k) 

fc '^'^kfj. 



where x = J2k<k^ ^k- Minimizing over x, 

,2 



Wo 



> 



> 



V- 2[v{0)+mmk<2k,v{k)] 2^0)' 



(3.17) 



(3.18) 



obtained by inserting 



v{0) + minfe<2fc^ v{k) 

in (3.17). Comparison with (3.13) shows that the energy 

density of any vacuum state is higher than the energy 
density of the ground state of the full diagonal model 
and, therefore, of the energy density of the QPGS. 

Proposition 2. Let v > 0. Suppose that the QPGS is 
unique. Then in the QPGS only mo can differ from zero. 

Proof. Uniqueness implies m_k = mk for all k. Sup- 
pose m_k = mk = n > for some k 7^ 0. Compare 
the energy i^m_k=mk=n of this state with the energy 
-E'm_k=o,mk=2n of thc State obtained by particle trans- 
fer from — k to k while keeping m_k' = mk' unchanged 
for k' ^ ±k. To compute the difference of these energies 
use 

Wk'-kmw = ^ «;_k',-km_k' = ^ Wk',kmk' 

k'T^itk k'ji±k k'T^itk 

(3.19) 

and 

v(2'k] 1 

w'k^-k - i«kk = 477^ + TjiMO) + v{2k)]xl > 0. (3.20) 



With them. 



2V V 



Em-^=Tnk=n ,,=i 



'm_k=0,mk=2ri 



k,k' 



= 2(«;k,-k - t«kk)n2 > ^^^n^ > (3.21) 

Thus the starting state is either not the QPGS or it is one 
but not the only one of them, contradicting uniqueness. 

Note that not only Wk,±k are positive but in fact all 
w^k,k' > 0. From Eq. (3.7) one easily deduces 

u(0)+u(k-k') ^;(0),^ ^ „^ , x „ 
Wk,k> > ^ ' '- + (/ik + /ik' + 2/ik/ik' ) > . 

(3.22) 

For (3.20) and (3.22) to hold, only i> > is necessary. 
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Theorem 1. Let > and w > 0. Consider the ba- 
sis (2.22) with Zk = 0. For /x > the QPGS is either 
the ground state of the full diagonal model or a small 
perturbation of it having the same energy density, or an 
eigenstate of Hqp with 



Let 



Xo oc Vv and mo oc W as V ^ oo, 

Xk - 0(1) for k 7^ , and xoV2 < (3.23) 

and with an energy density smaller than —^'^V/2v{0). If 
the decay of v{]s.) is slow enough {ip is sharp enough) then 
for small fi the last option holds. 

Example. If 

V?(r) = (47rA2)-it;(0) e""'/^^', t;(k) = u(0)e-(^'=)' 

(3.24) 

then the condition of sharpness is A < 0^/363 where 
as = v{0)m/ATrh'^, the scattering length in Born approx- 
imation. 

Remarks, (i) Because the average of A^o in a pairing 
product state {zq = 0) is 



(A^o) = ho + mo + 2homo, 



(3.25) 



(No) cx F in all cases. Thus, in the QPGS there is BEG 
of physical particles. 

(ii) When (3.23) holds for the QPGS, there is no macro- 
scopic and no generalized condensation for quasiparticles: 
{Mo) = mo cx VV and (Mk) = mk = for k ^ 0. More- 
over, the interaction being nonlocal, Swieca's proof [13] 
does not apply either. Therefore there is no reason to 
expect that the gap above the ground state of Hqp dis- 
appears. The gap will be derived in the next section. 

(iii) Of the two conditions i; > is the more important. 
> serves only to guarantee that the QPGS is not a 

vaciium state. It seems to be difficult to find a super- 
stable interaction whose QPGS in the VB scheme is a 
vacuum state. 

Proof of the theorem. From the two propositions above, 
the quasiparticle ground state energy is of the form 



Emo =wo + eomo -|- woo(toq - mo). 



(3.26) 



Because the last term is nonnegative, this can be smaller 
than the vacuum energy wq only if eo is negative. The 
minimum of (3.26) is attained with 



Too 



eo , 1 

- = my 



2woo 2 

that we suppose to be an integer, and its value is 



(3.27) 



Emo=mv =Wo- {woo - Bof = Wq - WoorUy (3.28) 

4W00 

We minimize Emo=mv with respect to xo which, in 
contrast to go, can take on any real value. For further 
use we introduce 



D 



V 



(3.30) 



Here and throughout the rest of the paper i/i and 1^2 are 
given by Eq. (2.8) with {Nh) = /ik and (^k) = Xk: 



(3.31) 



After some algebra one finds 



„2\2 r 







= {D + l-Cxo)il-x'oy[{D+l)xo + 2C] 

+1^0 (^1 + 1x1^ -(3.32) 

Since the second term is of order 1 while the first is 
typically of order V"^, the zeros of Gy are in a dis- 
tance 0{V~'^) to those of the first term. The first fac- 
tor, D + 1 — Cxq must be positive, otherwise eo > 0. 
Thus, Gv has three zeros in distance 0{V~'^) to ±1 and 
— 2C/(£'-|-l), respectively, and only in one of these points 
is the energy minimal. To minimize Em„=mv one has to 
maximize my, implying 1^2X0 < 0, see Eq. (3.34) below. 
For example, if 1^2 < then xo > 0, so the root we look for 
is either close to \2C/{D+1)\ or to 1. Now |2C/(L'+ 1)| 
may or may not be smaller than 1 and when smaller, it is 
the right choice. To see this note that if Xk = for k 7^ 
then C = and Gy (0) = 0. Thus, at xo = the energy is 
stationary and in fact minimum with value Emf- Indeed, 
(H^o) = in this case, therefore (H) = (Ho) > -Emf, see 
(2.7), (2.14), (2.17), and the lower bound is attained with 
the full diagonal model of ground state energy Emf- If 
Xk are nonvanishing and \2C / {D + 1)| < 1 then by con- 
tinuity it is still Xo = —2C/{D -|- 1) -|- o(l) that defines 
the minimum, and this remains true (with a negative o(l) 
correction) also if |2C/(i?-|- 1)| converges to 1 from above 
as V tends to infinity. 

Lemma 1. Under the condition that xq = —2C/{D + 
l) + o(l). 



lim 



1 



inf E„ 



2w(0)' 



(3.33) 



the energy density in the ground state of the full diagonal 
model. 

Proof. In terms of C, D and xq 



TOy 



(3.34) 



Substituting the assumed value of Xq and ^/l — Xq = 
1/(1 + 2ho), 

{No)=ho + {l + 2ho)mv = ^^{n-ui) + o{V). (3.35) 



8 



Adding to this the number of uncondensed particles, 



(AT) 



v{0) v{0) 



L^^vO^)hu + o{V). (3.36) 



Suppose now that is the energy- minimizing set. Then 



dfj. 



-{N) 



±^^Y.<^)K + o{V) (3.37) 



where the o{V) correction can be negative. On the 
other hand, at /z = we have 5k = 
and Emo=mv{{9Vs) — (because y > 0). Integrating 
Eq. (3.37), 



En 



'2v{Q) 



(3.38) 



Because this is the energy minimum, 

£^mo=niv ({.9k}) ^ ^-MF- 

The second term in the right member of Eq. (3.38) is 
nonnegative and, thus, of oiV) [implying also ^k = o(l)], 
otherwise the above inequality could not hold. This, how- 
ever, proves Eq. (3.33). 

An energy density smaller than —fj,^ /2v{0) can only be 
obtained if the energy-minimizing (/k for k ^ are such 
that \2C/{D + 1)1 > 1 and remains separated from 1 as 
V tends to infinity. In what follows we shall investigate 
this possibility. 

Lemma 2. Suppose that |2C/ (Z? + 1) | > 1 and remains 
separated from 1 as V tends to infinity. Assume also 
/Lt — z^i + l^al > 0. Then the energy- minimizing xo and 
also my are of order ^/V. Furthermore, 



inf E„ 

go 



V 



3t;(0) 

where J-,o is wq at = 0, 



{^l-Vl+\u2\f+J^o+o{V) (3.39) 



2V 



lk#0 



+ W E ^^(k-k')(XkXk'+/ik/ik') • (3.40) 

k.k'^iO 

Proof. Under the condition of the lemma and with the 
choice C < 0, xo is necessarily in o(l) distance to 1. To 
compute its asymptotic form, in Gy we insert xq = 1 
everywhere except for I — Xq and obtain to leading order 

Gv{xo) = g-(D+l+|C|)(2|C|-D-l)(l-x2)2. (3.41) 



Solving Gv{x(j) = yields 

V(i^ + i + |c|)(2|c;|-D-i) 

This shows that a;o = 1 - 0(1^"^), 

■{D+l + \CmC\-D-l) 



(3.42) 



Xo 



27 



oc VV (3.43) 



and with xo = xa/{2y^l — Xq) 



mv 



1 D + 1 + 
"2^ 3X0 



\C\ 



oc 



(3.44) 



From here the form (3.39) of the energy is obtained by 
substitution into Eq. (3.28). 

To complete the proof of the theorem it remains to 
show that (3.39) can be smaller than Emf and the dif- 
ference can be of the order of the volume. Let 

ui = i/o + v{0)p', uo = ^^v{\i)hk, p'=^Yhk 



V . 



and 



k^^O 

(3.45) 



1 

V' 



J^O = JjJ^O+l^p' = :^^£(k)/lk + ^{p'f 



'(0), 



^ J2 t'(k-k')(xkXk'+/ik/ik') -(3.46) 



21/2 



Now 



k,k'#0 



inf ErnQ=mv 

V go 



2v{i)) 



0(1) 



1 



{li''-2[v{0)p' + 2{\u2\-vo)]iJi 



+6v(0)i^o - 2(|i/2| - vif] 
is negative if 

A = [v{())p' + 2{\u2\ - vo)f + 2{\v2\ - vxf - 6i;(0)j^o > 
and 



(3.47) 



(3.48) 



u(0)p'-h2(|z/2|-z^o)-VA <p.< t;(0)p'-h2(|zy2|-!^o)+VA- 

(3.49) 

The conditions eo < (or my > 0) and 2\C\/{D + 1) > 1 
of the validity of (3.39) imply 

fi -W2\ = 

v{0)p' + 1^0- \v2\ <p< v{0)p' + 2\u2\ + yo 

= vi + 2\v2\- (3.50) 

The intervals (3.49) and (3.50) for p clearly overlap if 
\u2\ > fo which can be realized if all Xk have the same 

(negative) sign because /ik < |xk|- We therefore do not 
consider the otherwise interesting possibility that Xk can 
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vary in sign. On the other hand, for interactions that are 
bounded at the origin Xk < for all k 7^ implies that 
A becomes negative for large /x. Indeed, with Eqs. (3.45) 
and (3.46) the discriminant can be rewritten in a more 
convenient form, 



A = 6(N-.o)^-M^ 



^£(k)/lk 



^ t;(k-k')(xkXk'+/ik/ik') (3.51) 



k.k'T^O 

Using the Schwarz inequalities 

^ Y: .(k-k')xkXk'>-; 



y2 



^ V{k - k')/lk/lk' > 1^0 



(3.52) 



k,k'7^0 



one obtains 



A < 3(|zv2| - i^of - 61^0^21 - ^ ^£(k)/ik. (3.53) 



Now 



kjtO 



because /ik < |Xk| < /ik + ^, and therefore 

A < ^<p{Of - 61/0^2! - ^ (3.55) 



It is then clear that the positivity of A places an upper 
bound on |xk| and, via (3.49) or (3.50), on jj,. 

We now show that under some additional condition 
on the interaction, for any small enough /i the difference 
(3.47) can be made negative. Because (p > and is in- 
tegrable, w(k) is continuous and |i'(k)| < v{0). Choose 
< 7 < 1 and fc^ > such that v{0) > v{k) > 7t;(0) for 
k < kj. Choose moreover Xk < and 



h ifO <k< k^ 
if fc > A;-v 



(3.56) 



Let 



^7 = ^Ei-^- 



Then with x = y^h{h+l), 

1 . 3 
- ^£(k)/ik « -/iB^e(k^) = 



ahk^ 
IOtt^ 



(3.58) 



{a = fr/2m) and 



we obtain 



The lower bound is positive for small enough h if 

5(27^ - l)v{Q)k^ > (67r)V. (3.61) 

This shows that 7 > 1/^/2 must be. For example, with 
the choice 7 = 0.9 and k^ = {A-k)^ a / v{Q) the condition 
on the interaction reads 



7(k) > 0.9t;(0) for k < 



(47r)V 
^;(0) 



(3.62) 



The positivity of the r.h.s. of (3.60) imposes x < l/vl2, 

but no lower bound on x- Therefore the (positive) lower 
edge of the interval (3.49) can be arbitrarily close to zero. 
Similar holds to the lower edge of the interval (3.50) if it 
turns out to be positive. Thus, under condition (3.62) the 
difference (3.47) is indeed negative for any small positive 
II. Although ^5(0) < OC' for the Gaussian example^ (3.24), 
Eq. (3.62) is compatible with the non-integrability of v 
and, hence, with (p{0) = 00. 



B. Gap to excitations in the VB scheme 

Adding n quasiparticles of wave vector k = to the 
QPGS needs an energy 



E„ 



while the gap to k 7^ excitations is 



(3.63) 



E 1 



^mo—mv 



WOk 



Ck + 2wokTOy = ek ^eo + wok- (3.64) 

Woo 

Equations (3.14)-(3.15) are special cases of Eqs. (3.63)- 
(3.64), obtained with all gk = 0. If the QPGS is charac- 
terized by Lemma 2, the gap to k = excitations is 



Woo 



KO) 
2V 



(1 + 6x0) 



MO) 



-Xo 



-(/i + \u2\){-n + + 2\u2\). (3.65) 



(3.57) To compute (3.64), E,, 



has to be minimized w.r.t. 
Xk- The analogous procedure in the NStJAW scheme 
will be discussed in detail in the next section. Here we 
present only the limit of (3.64) as k — > [19]. It can be 
obtained from 



2{H-Vi)-\U2\ 

lim Xk = — |- 

k^o n — Vi + 4\t/2\ 



(3.66) 
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valid for the miniinizing a;k, cf. Eq. (3.29). The result is 

lim (-£/77io=mv)7^k=l -^mo=mv) 
k— >0 



= Y - + + 1^1+ 5\U2\). (3.67) 

This is somewhat higher than wqq, the gap to k = 
excitations. The discontinuity of the gap comes from the 
discontinuity of a;k at k = 0. Indeed, from the strict 
inequalities (3.50) one obtains 

lim Xk<xo = l- 0(y"M. 



IV. THE NOZIERES-SAINT 
JAMES-ARAKI- WOODS SCHEME 

A. QuEisipEirticle ground state 

Comparing Eq. (3.39) with the lower bound (2.14) one 
observes a prefactor | in the former instead of ^ in the 
latter. The source of the difference is D'^{Nq) which was 
replaced by when (2.14) was derived but which is 



D\No)^-{Nof^2xWy 



2V^ 
9v{0y 



-{iJ. -Ui + \l/2\f 

(4.1) 

in case when (3.39) is valid. In the VB scheme the fluc- 
tuation of the number of condensed particles is anoma- 
lously large, manifesting itself in an excess energy den- 
sity. In 1982 Nozieres and Saint James [10] studied a 
modification of the VB scheme in which they dropped 
the Bogoliubov transformation in the zero mode and 
applied instead Bogoliubov's c-number substitution for 
oq. In the present work this corresponds to setting 
go = and applying a shift (2.18) only in the zero mode: 
«o = 0^0 ~ = ^Oj Ok = Ik for k 7^ 0. Because then 



D^No) = (No) = \z\ 



(4.2) 



by minimizing the energy w.r.t. z, up to a correction of 
order 1 the lower bound (2.14) could be attained. As 
Eq. (2.15) suggests, the minimizer was 



V 



(4.3) 



see below. Thus, remaining within the family of pairing 
product states, only an improvement of order 1 can be 
gained on the NStJ result. It is nonetheless instructive to 
see the origin and some consequence of this improvement. 
Before arriving there, first we will have a closer look at 
the NStJ setting. Nozieres and Saint James considered 
only the vacuum state $o. Here we take into account all 
states of the basis (2.22) and prove the following. 

Theorem 2. Let (p > and v > 0. Consider the basis 
(2.22) with go = and Zk = zS]^^. Under the condi- 
tion of uniqueness, the QPGS is a quasiparticle vacuum 



state with an energy density smaller than — /i^/2?;(0) for 
all /X > 0. If (^(0) = oo, for /i large the full parti- 
cle density is p = [n + o(/x)]/w(0) where the o(/x) term 
tends to infinity with /z. If ip{0) < oo, for /x large 
p = + (p{0)/2]/v{0) + o{iJ,°). In all cases the uncon- 
densed density tends to infinity with jj, but slower than 
/U, e.g. p' ~ fj,^^^ if the decay of v is fast enough. 

Remarks, (i) In the ground state of a system of clas- 
sical particles whose interaction satisfies w > 0, f (0) > 
and z;(k) = for fc > /cq, the exact relation p = 
[jj, + (p{0)/2]/v{0) holds above a threshold density ~ fcg, 
see Ref. [20]. 

(ii) The role of (/J > is to guarantee that the QPGS is 
a vacuum state - just the opposite of its role in the VB 
scheme. 

Proof. The argument given in Proposition 2 remaining 
valid, uniqueness implies that in the QPGS only toq can 
be nonzero. Let denote such an element of the basis 
and {B) the expectation value of the operator B in this 
state. For k 7^ from (3.8)-(3.10) we find (TVk) = /ik, 
(^k) = Xk and D'^{Nu) = xl- Because = 0, {H^o) = 
J-,0, cf. Eq. (3.40). For k = now we have 



Too, 



(4.4) 



and 



D^{No) = \z\''{l + 2too) = {No) + mo{2\zf - 1). (4.5) 

Choose z real and 1^2 < 0- Insert {No)-{Ao) and D^{No) 
from above into (2.9) to find 



{H) = J^o-{^i-i^i + W2\){No) 
+mo ( 11^2 



^-^{Nof 
2V 

vi0)^2 v{0) 



V 



2V 



■ (4.6) 



If \u2\ is of order 1, the coefficient of toq is positive for 

any z, and the minimum of (H) is obtained with Too = 
and, provided that 11 — i^'i + \i'2\ > 0, with (No) = z^ 
given by (4.3). Its value is 



Eq = min($05 H^o) 

z 

V 



V 



2v{0) 



{p. - lyi + \p2\f + J^o 



2w(0) ~ ~ ^2)^ - (z^o + i^2)p'V + ^ £(k)/ik 
+W E ^'(k-k')(XkXk'+/ik/^k'), (4.7) 



2V 

k,kVO 

see Eqs. (3.45), (3.46). Note that 



p={N)/V = po + p' = 



^2 



(4.8) 



thus, the first term of the last member of (4.7) is 
—v{Q)p'^V/2. The minimum of £^0 over {fifklkT^o is to 
be compared with the lowest energy one can obtain if 
the coefficient of toq in (4.6) is negative. 
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We show that if the coefficient of mg in (4.6) is neg- 
ative then the energy is at the best by order 1 below 



and 



-/Lt^V/2w(0). In this case < \ and \i'2\ < 



(The 



latter may occiu" cither because Xk fi^rc small or because 
they oscillate in sign, in which case v\ may still be of or- 
der 1.) Because of the superstability of the interaction, 
in the QPGS mo can be at most of order V . Therefore, 
the contribution to (i?) of the term proportional to mo is 
at most of order 1. For the same reason |z^2|(-^o) — 0(1). 
As for the rest, there are two cases to be distinguished, 
(i) \i — v\ < 0. Then the minimum of (7J) is, apart 
from a correction of order 1, that of J^q which cannot be 
smaller than -ii^V jlvifS), cf. Eq. (3.16). (ii) /i- i^i > 0. 
Minimizing the right member of Eq. (4.6) with respect 
to (iVo), 

(if)+0(l)> J.o-^(M-'^i)'- (4.9) 

Substituting J-,o from (3.40), j^i = i^o + f (O)p', and ap- 
plying the Schwarz inequalities (3.52), after simplification 
we find 

(if) + 0(i) > ^£(k)/ik + :jJ:T(M-^(0)p')^o '''^ 



«(0) 



> - 



2w(0) 



2t;(0) 
(4.10) 



because \x — v{0)p' > uo> 0. 

On the other hand, by choosing Xk < and |z/2| small 
but non- vanishing as V tends to infinity, for any > the 
energy (4.7) can be made smaller than — /i^V/2t;(0) by an 
amount of order V. Indeed, in this case —i/q — ;/2 > 0, cf. 
Eq. (3.54), so an improvement comes from the first term 
of (4.7); the last three terms are positive, independent 
of fi and can be made arbitrarily small by decreasing 
IXkl- Consider, for example, the choice (3.56). Because 
ft < X^) in the limit when x ^ we can drop terms of 
order or smaller in (4.7). Therefore 



-7B^XM+0(x' 



V 2vi0) v{0) ^ > - 2v(0) 

(4.11) 

While (4.11) is the relevant formula for the energy den- 
sity if /i is small, the situation is quite different in the 
other asymptotic region, when /x tends to infinity. It 
is important to realize that Xk minimizing the energy 
do not remain small when /z increases, as the argument 
leading to (4.11) would suggest it. On the contrary, they 
tend to infinity with /U, although more slowly. This will 
become obvious when we minimize (4.7) in the variables 
Xk- Prom Eq. (4.7), 



dEo ^ 
where 



2a;k [v (k)po + I'ok - I'd - + £ (k)] -|- 2 [i> (k)/9o + i^2k] 

(4.12) 



Po = 



(No) 11-V1-V2 



V 



v{Q) 



(4.13) 



^o\^ = h^ ^^(k - k')/ik', z^2k = ^ u(k - k')xk' 

k'#0 k'#0 

(4.14) 

(i^oo = z^o, 2^20 = z^2, fi-Vi = i^ik-) Then dEo/dxk. = 
yields 



v{k)po + 



v{\s)pq + I/Qk - I'D - f2 + £(k) ' 



(4.15) 



This is still an implicit relation because po,yo-, i^ok) ^2-, t^2k 

depend on all Xk'- Let /i tend to infinity and assume 
that Ixk'l do not increase. Because po tends to infinity, 
^ —1, i.e., Xk — * —00, contradicting the assumption. 
So Ixk'l do increase and therefore the imcondcnscd den- 
sity p' and all the v diverge with a diverging p. To allow 
a;k to tend to —1, |fjk|/po must tend to zero. Moreover, 
because fok and p' are of the same order, p' /po tends to 
zero as p. tends to infinity. Comparison with (4.13) then 
shows that po is of the order of p and p' = o{p). The 
result for p follows from Eqs. (4.8) and (3.54) with the 
remark that 



(4.16) 



as /ik ^ oo- 

More precise asymptotics in p can be obtained if we 
choose /ik ~ \h{k) and Xk < where ft.(k) > satisfies 
/ £(k)/i(k) dk < 00 and / v{k)/h{k)'^ dk < 00. For such 
an /i(k) to exist v must decay fast enough; for example, 
a Gaussian interaction is suitable. Now A is a parameter 
to minimize Eq. As /U — > 00, the minimizing A tends 
to infinity as well. The terms of Eq/V determining the 
asymptotic form of X{p) are 

pj[v{k)/h{k)]dk A2 



(47r3)?;(0)A 



(2 



J J w(k-k')/i(k)/i(k')dkdk' 



where the first term comes from the expansion (4.16). 
The minimum is attained at 



A: 



pj[vik)/hik}]dk 



2v{0) JJv{k- k')/i(k)/i(k') dkdk' 



(4.17) 



We find therefore p' = (A/Stt^) / h(k)dk ~ pi. This 
finishes the proof of the theorem. 



QuasipEirticle Hamiltonian and gap to 
excitations in the NStJAW scheme 



With the help of the basis (2.22) one can again define 
a quasiparticle Hamiltonian. The contribution of iJ-,o 
to Hqp is the same as in the VB scheme, the difl^erence 
comes from a difi'erent contribution of Hq. Compared 
with (3.3) the changes are as follows. First set ho = 
Xo = at every occurrence. This yields 



U!q = J^o, eo = -p + iyi, Woo 



<0) 
2V ' 



(4.18) 
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ek = (l + 2/ik)(£(k)-/i + z/ik) + 2xkZ^2k (4.19) Then 



and 



wok = ^{l + 2hi,)[v{0)+vik)] (k^O) (4.20) 



where 



i^ik = v{0)p' + z/Qk 



(4.21) 



{vio = v\, z^i,_k = z^ik). Second, replace wq, ek^^o and eo 
respectively by 

iZJo = «;o - (m - + (Rez2)z^2 + ^|^|^ (4.22) 



ek = ek + ^(1 + 2ftk)[t;(0) + t;(k)] + %i;(k)R, 



and 



eo = eo + -|^|2:r- 



(4.23) 



(4.24) 



The other cocfScicnts arc unchanged. Because the 
ground state of ifqp is now $0) adding a quasiparticle 
of wave vector k to the ground state costs an energy 
Ck. Inserting z real from Eq. (4.3), using Eq. (4.13) and 

Ck = (l+2ftk)[t;(k)po(l+a;k)+z^ok-z^o+a;kZ^2k-z^2+e(k)] 

(4.25) 

and 



eo = /i - 1^1 - 2i/2 = ■i;(0)po + k2|. 
Prom (4.15) 

w(0)po +1^2 M - i'l 



lim Xk = — - 



(4.26) 



(4.27) 



k^O i'(0)po — Vl II — V\— 2V2 

see the remark [19]. Substituting this into ek and using 
1 + 2/ik = l/Vl - a;k we find 



^im ek = 2^v{0)po\v2\ > 0, 



(4.28) 



the same as the gap vahie in Refs. [15] and [10]. As 
expected, limk^o Ck 7^ eo because limk^o a;k 7^ a;o = 0. 
In effect. 



< lim ek < eo. 

k^O 



(4.29) 



The discontinuity of the gap at k = will be removed in 
the next section. 

In terms of the numerator and the denumerator of a;k, 
eik has a nice form. Let 

ak = v{}s)P0 + Z^2k, /3k = w(k)/9o + Z^Ok - Vq - V2 + £{\^) ■ 

(4.30) 



a;k = -ak//3k, ek = V(/3k + Q!k)(/3k - ^k), (4.31) 
Choose now k 7^ such that 

v{Q)p' «: e(k) «: v{Q)pq. (4.32) 



Keeping only the largest term, 

ek « V2^^(k)po£(k). 



(4.33) 



This is Bogoliubov's dispersion relation that was repro- 
duced ahcady by Girardeau and Arnowitt [15] and is 
obtained lierc within the NStJAW scheme. With energy- 
minimizing Xk the interval (4.32) does not extend down 
to zero because the optimal p' is separated from zero. 
Hence, the linear dispersion exists only as an extrapo- 
lation near k = 0. This, of course, is consistent with 
Eq. (4.28). As p diverges, the interval [v{Q)p' ,v{{))po] 
widens and its lower bound shifts upwards e.g. as /u^/^, 
cf. Theorem 2. As p tends to zero, the interval gets closer 
to zero. Because of the relation v{Qi)pq = p — vi + \v2\, 
the decrease of p is consistent also with po fixed and w(0) 
and all ft-k, |Xk| and, therefore, p' decreasing. It is in this 
sense that Bogoliubov's result for a weakly interacting 
Bose gas with an almost complete BEG can be recov- 
ered. When we take the limit w(0) 0, Eq. (4.15) shows 
that Xk ^ (and, thus, (iVk) 0) for any fixed nonzero 
k. However, if k decreases with w(0) so that e(k) remains 
in the interval (4.32), then Xk does not tend to zero; ac- 
tually lim|xk| > 1/2 in this joint limit. Note that the 
joint limit can be realized only in infinite volume, when 
k can vary continuously, and the result does not contra- 
dict p' — > 0. In finite volumes, for u(0)po < I? j^ml? 
there is no e(k) in the interval (4.32), and Xk eventually 
tends to zero with v{^) for every nonzero k. 



C. «o 7^ and ffo 7^ 

Here we discuss a slight improvement of the NStJAW 
scheme. Comparison of Eq. (4.7) with the lower bound 
(2.14) shows that by letting 50 be nonzero and by mini- 
mizing the energy with respect to it, only a gain of the 
order of unity can be attained. Minimization in 50 is, 
nonetheless, interesting because it makes the excitation 
energy continuous at k = 0. 

Choosing nonzero affects only (iJo) which can be 
computed from Eq. (2.9). The relevant values to be sub- 
stituted here are 



(TVo) = + (iVo), (^0) = + (^0) (4.34) 



and 



i?2(iVo) = z2[l + 2(iVo)] + 2z2(Ao) (4.35) 
where z is real and 



(iVo) 



(ao«o, 



= ft-o + (1 + 2ho)mo, 

Xo + 2xomo. (4.36) 
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These are averages in a state (2.22) with Zk = = 

for k 7^ and Zq = z; ho and xo are related to via 

Eq. (3.2). Substituting (4.34) and (4.35) into Eq. (2.9) 
we find 



(Ho) 



11-VI-U2 
+ 



''^^\{No) + {Ao)]] {No) 



V 



2V 
2V 



{Nor - ^2 [{No) - {Ao)] 
(iVo)[l + 2(7Vo) + 2(io)] 



(4.37) 



The minimum of {Hq) as a function of (A^o) reads 
v{0) 



{Ho)=- 

obtained with 
{No) = Vpo = 



2V 
2V 

V 



{No)^ - U2[{No) - {Ao)] 
(iVo)[l + 2(iVo)+2(Io)] (4.38) 



{n-u,-i^2)-{No)-{Ao). (4.39) 



The TOo-dependent part of {Ho) in (4.38) reads 

v{0) 



(1 + 2ho) 



\U2\{1 - Xo) 



2V 
V 



(l + 4/io)(l + a;o) 



mo 



(l + 2/io)^(l + a;o)m^. 



(4.40) 

The difference Eo — {Hq), where Eo and {Hq) are taken 
from Eqs. (4.7) and (4.38), respectively, cannot be larger 
than of order V'^. This, however, means that {No) and ho 
and mo together with it are of order , implying that 
the expression (4.40) is nonnegative for V large enough 
and its minimum is attained at mo = 0. Setting mo = 
we obtain the vacuum expectation value 

^0 = ($0, H^o) = J^o ^ 



2v(0) 



+(/" - i^i)(Xo + ho) - 2u2ho 



V 



ho{l + 2ho + 2xo). 



(4.41) 



Here the second line contains the correction to (4.7) 
which is negative if xo < 0. For minimization consider 



dEo 



Xo 



jj.- Ul - 2V2 



V 



(1 + 4/io + 2xo) 



+ 



dEo 1 9x0 = yields 



xo = - 



^ho 
1/ "-0 



2v{0) 
V 



ho 



(4.42) 



^-ui-2u2-^{l + 4.ho + 2xo) 
viO)po + V2~ ^{hg - Xo) 
t^(0)po-i^2- 4^(1 + 3/10+ xo)' 



(4.43) 



In computing dEo/dxv for k 7^ 0, (4.12) must be re- 
placed by the derivative of (4.41). We find 



9xk 9xk 



90=0 



dEo dEo t;(0) + Kk) 4t;(k) 

2a;k y (Xo + ho) y~ho 

{AAA) 

where the first term is given by (4.12). Equating this 
with zero we obtain 



v{\£)po + I/2k - ^(/io - Xo) 



w(k)po + i/Qk -V0-V2+ e(k) -^{ho + Xo) 



(4.45) 

Note that in Eqs. (4.43) and (4.45) po is substituted from 
(4.39). It is seen that 



lim Xk =xo + 0(V"^), 

k— >0 



(4.46) 



cf. remark [19]. 

Consider now the associated quasiparticle theory. Hqp 
is obtained from (3.3) by replacing the VB coefficients wo, 
Bkjio and eo (all defined with go ^ 0) with the expressions 



Wo 



wo-{p-yi- V2)z' + + --^7?{2ho + xo) 



.2 , 'K0).4 , ^).2, 

V 

z\. . „ 2xk 



2V 



ek = ek + y(l + 2/ik)b(0)+t;(k)] + =^t;(k)^2 



eo = eo + ^^z^(l + 2/io + Xo) 



(4.47) 



where z is real, and keeping the other coefficients un- 
changed. Now (A^o) = + ho, therefore z^/V = 
Po + 0(y-^). Thus, with (4.39) 

ek = (l-|-2/ik)[v(k)/)o(l + a;k) + i'ok-fo 

+XkU2k - U2 + s{]f.)] + 0{V-^) (4.48) 

while 

eo = (1 + 2ho)[viO)po{l + xo) - - xo)] + 0{V-^). 

(4.49) 

Comparing the above two equations we find 



lim Ck = eo+0(F-i) = 2^v{0)po\u2\+0{V-^), (4.50) 

k— >0 

the continuity of ek at k = 0, up to an error of order 



V. SUMMARY 

In this paper we have investigated variational wave 
functions for the ground state of interacting homogenous 
Bose systems. We have not invented new trial functions 
but studied with much care those introduced by Valatin 
and Butler [9] and modified later by Nozieres and Saint 
James [10] — though we completed the trial states of the 
latter into a basis in Pock space. Both sets of trial func- 
tions are quasiparticle product states grouped in bases 



14 



in Fock space, each basis corresponding to a choice of 
one parameter for every pair {k, — k}. Their difference 
is in the choice of this parameter for k = 0, associated 
with Bogoliubov's canonical transformation for VB and 
with a coherent state for NStJ. Each basis consists of a 
quasiparticle vacuum state and excited states with an ar- 
bitrary finite number of quasiparticles. However, it is not 
a priori guarantied that the vacuum state has the low- 
est energy within a basis, therefore energy minimization 
must extend to the full basis for a given parameter set. 
To solve this problem, one has to set precise conditions on 
the interaction, which in our case are positivity of both 
the pair interaction and its Fourier transform. Assuming 
that the overall energy minimum is attained on a sin- 
gle state, we have shown that this state could contain at 
most quasiparticles of zero momentum. It then turns out 
that in the VB scheme the lowest energy state is never 
a vacuum state, as was supposed by Valatin and Butler. 
For slowly decaying w(k) and small it is a state with 
~ VV quasiparticles, and we could not exclude the possi- 
bility that in other cases it is an eigenstate of the number 
operator No, namely, the ground state of the correspond- 
ing mean field model defined by WMF(k) = w(0)(5k.o- In 
the first case both the mean and the standard deviation 
of the number of physical particles of zero momentum 
are proportional to the volume. This anomalous fluctua- 



tion causes an excess energy density, not present when a 
coherent state is used in the zero mode. Therefore, the 
lowest attainable energy density in the NStJAW scheme 
is always below the minimum obtained in the VB scheme, 
and is provided by a vacuum state. 

Unbiased minimization of the energy leads to a gap 

in the spectrum of the associated quasiparticle Hamilto- 
nian. We have computed these gaps both at k = and 
in the limit k ^ 0. The two values do not coincide in 
either schemes. In the case of the NStJAW scheme we 
have shown that the continuity of the gap at k = can 
be restored through an additional energy-minimizing Bo- 
goliubov transformation in the zero mode. We have also 
discussed how Bogoliubov's linear dispersion relation for 
a weakly interacting and almost completely condensed 
Bose gas can be obtained within the NStJAW approach, 
in spite of the gap in the spectrum. 
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